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SYMBOLS 


EQUATION OR 

BASIC QUANTITIES AND COEFFICIENTS FIGURE 

A The dimensionless coefficient indicating the factor by 

which (5p/p) is multiplied to yield (5T/T) (20) 

G The standard geopotential gravitational constant 

9.80665 m 2 s' 2 (mV (1) 

H Scale height in units of geopotential height, m' (or km') (1) 

h Geopotential height in units of m' (nearly equal to geo- 

metric height, but accounting for height variation of the 
acceleration of gravity; that is, Gdh = g(z)dz) (1) 

k A dimensional constant 0.0341632 K/m' resulting from 

the combining of the three constants GM/R into a single 
value (4) 

M Mean molecular weight of air 28.9644 kg kmol' 1 (1) 

p Atmospheric pressure N m" 2 (1) 

R Universal gas constant 8.31432 X 10 3 joules K' 1 kmol" 1 (1) 

R Ratio of S _ to S Figure 6 

q.max q,min 

S The value of a particular series of terms involving 

density ratios (30) 

T Atmospheric temperature K (1) 

x A general designation for a function of a number of 

generalized variables ( 1 4) 

y A general designation for the set of variables of which 

x is a function (14) 

a A subscripted dimensionless coefficient with a value 

near unity (23) 

p Atmospheric density kg m" 3 (7) 

iii 



mil 


EQUATION OR 

SINGLE SUBSCRIPTS FIGURE 


1 

Refers to a particular height h x so that p x , T x , p x ,and h x 
in a single equation all refer to conditions at the same height 

(4) 

2 

Similar to 1 

(4) 

i 

Refers to an individual member of a set of related 
variables having the general designation y. 

(14) 

i 

Also refers to a variable in an isothermal atmosphere as p. in 
the graph of (p./p r ) versus h in figure 1 

(14) 

j 

A general designation for an integer which may vary be- 
tween 2 and q, and which is simultaneously associated with 
a geo potential height as lv and with the related density as p j 

(21) 

q 

A general designation for an integer which may have any 
positive value, and which is associated with the lowest 
density-height data point (that is, h q , p q ) involved in 
a particular evaluation of an integral to determine the 
value of the related temperature T q 

(11) 

q 

Also used as S q where q implies the number of terms 
in the series S 

(30) 

r 

Designates a specific reference value for the basic 
quantity of which it is a subscript, that is, p r and p r 

Figure 1 

DOUBLE SUBSCRIPTS 


1,2 

Refers to a particular height layer as between h 2 and h x , 
where h x and h 2 are ordered in such a way that the 
quantities within the subscripted parentheses have par- 
ticular signed significance; that is, Ah is always positive 
while AH is positive only when dT/dh is positive in the 
related layer 

(4) 

H,j 

Similar to 1 , 2 

(25) 

jj+i 

Similar to 1 , 2 

(26) 

q-l,q 

Similar to 1 , 2 

(27) 


iv 


EQUATION OR 
FIGURE 


q, max Designates maximum value of S q Figure 6 

q, min Designates minimum value of S q (36) 

OPERATORS AND FUNCTIONS 


— 

Overbar as in T or H indicates the mean value of T or H 
for the related layer 

(4) 


Square root as^/T* 

(20) 

/ 

Integral as / p(h) dh 

(ID 

a 

Partial differential as in 9x 

(14) 

d 

Differential as in dh 

(1) 

£n 

Natural logarithm as Cnp and £np 

(1) 

A 

An increment as in Ah and AH 

(4) 

5 

An increment or random uncertainty as in 5y., Sp x , 
6Ah, 5Tj, 5p., 5T q , and 6p q 

(14) 

2 

A summation of a set of terms 

(12) 
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AN ANALYSIS OF THE ERRORS ASSOCIATED WITH 
THE DETERMINATION OF ATMOSPHERIC 
TEMPERATURE FROM ATMOSPHERIC 
PRESSURE AND DENSITY DATA 


R. A. Minzner 

Goddard Space Flight Center 


INTRODUCTION 

Relationships between pressure, temperature, and height in the earth’s atmosphere are well 
known and for many years have been the basis for height determination in balloon- 
radiosonde flights in which pressure and temperature were measured during meteorological 
probings of the lowest 30 km of the earth’s atmosphere. When rocket vehicles extended 
the potential atmospheric probing capabilities to heights up to 1 00 km, into regions where 
pressure is still measurable, but where existing technology does not allow for immersion sen- 
sing of temperature, these same mathematical relationships were used to extract temperature, 
from the measured pressures and radar-derived rocket-height data. The unsmoothed temper- 
ature-height values from many of these soundings represented a very jagged height profile, 
with the degree of jaggedness apparently increasing as the height increment between suc- 
cessive data points decreased. Investigators usually have been unable to determine how 
much of this jagged ness represents real temperature variability and how much is attributable 
to measurement error. 

The hydrostatic equation and the equation of state lead to another set of height relation- 
ships: those between atmospheric density, temperature, and height, such that temperature 
can also be computed from density-height data. In cases where such computations have 
been made, particularly in the height region of 30 to 100 km, the resulting temperature- 
height profiles appeared to be less jagged than those derived from pressure-height data 
with comparable height resolution. 

The apparent difference in the jaggedness of density-derived temperature-height profiles 
from those associated with pressure-height data suggest that the height increments of the 
pressure and density data do affect the uncertainty in the derived temperatures, and that 
the influence of height increments in relation to density data may be different from that 



in relation to pressure data. Obviously, error analyses, which involve both the pressure- 
temperature-height relationship in one instance and the density-temperature-height relation- 
ship in another instance, are needed. The error analyses presented in this paper confirm the 
fact that important differences do exist between these two sets of relationships, particularly 
in regard to the influence of the height increment on the propagation of measurement un- 
certainties into the temperature-height profile. These differences strongly favor the use of 
density-height data over pressure-height data. 


FUNDAMENTAL CONSIDERATIONS 


Pressure-Height Relationships 

The equation of state, when combined with the differential form of the hydrostatic equation 
to eliminate density p, yields an expression frequently referred to as the hypsometric 
equation: 

dEnp = -GM = -l_ 
dh RT ~H 

where 

p = atmospheric pressure, 

h = a measure of the height above sea level, in geopotential meters m' 

T = absolute temperature of the atmosphere at h, 

R = the universal gas constant, 

M = the mean molecular weight of air, 

G = a constant when h is expressed in geopotential, and 
H = the scale height in geopotential units. 

Solving equation (1) for T and H, respectively, yields 


and 


GM -dh 
R dfinp 


( 2 ) 


RT _ -dh 
GM d£np 


( 3 ) 


When values of £np versus h are known from numerical data rather than from an analytical 
function, it is convenient to replace the expression -dh/d£np, which applies to a specific 
height, with a numerical approximation (h 2 - h t l/lTnPj - £np 2 ). In this approximation, 

Pj is the pressure at hj , and p 2 is the pressure at h 2 . The approximation therefore represents 
the mean value of the reciprocal of the derivative over the height interval (h 2 - h t ). When 


2 



the point value of the derivative in both equations (2) and (3) is replaced by the numerical 
approximation, the related values of T _and H no longer apply to a single height, but 
rather become the mean values T and H, respectively, associated with the height interval 
(h 2 - hj ). Thus, we have 


T = 


GM 

R 


h„ 


= k 


( Ah )j ^ 


(4) 


£np t - Cnp 2 Cnpj - Cnp 2 

where k = GM/R = 0.0341632 K/m', and (Ah)j 2 = h 2 - h t , and where, from the relation- 
ship between T and H implicit in equation (1), we see that 


h 2 -hj (Ah) 1 2 

Knpj - Knp 2 Cnp t - £np 2 

From equation (5), we obtain the following expression which will be of special importance 
in the section on Uncertainty of Temperatures Deduced from Pressures: 

— 1 - - (!) 

Cnpj - Cnp 2 yAh/j 2 

The interrelationships between £np, dCnp/dh, T, and H, all as a function of height as ex- 
pressed by equations (1), (2), and (3), are shown pictorially by four of the six height-related 
graphs of figure 1. This figure, which depicts the properties of a portion of the U.S. 
Standard Atmosphere, 1975* was developed from a preliminary set of abbreviated tables 
representing a portion of this revised standard atmosphere (Kantor and Cole, 1973). Geo- 
potential height is scaled linearly along the abscissa of figure 1 . Atmospheric pressure, 
which is plotted in the form of the natural logarithm of the ratio (p/p r ), has its ordinate 
scale at the left of the figure and is depicted by the lower of the three lines diagonally 
crossing the entire figure from upper left to lower right. In this presentation of pressure, 
p r in the ratio (p/p r ) has been made equal to one newton per square meter (1 N m' 2 ) so 
that the numerical values of (p/p r ) are those of p in N m' 2 , while the scale of Cn(p/p r ) at 
the left of the figure implies no dimensions. 

The graph of dCnp/dh shows that the slope of 8n(p/p r ) versus h, or equivalently the slope 
of Cnp versus h, is constant only in height regions over which T is invariant. (The value of 
dfinp/dh is negative and inversely proportional to T.) These portions of the graph of 
dCnp/dh versus height corresponding to constant T (the three height intervals 1 1 to 20 km', 
47 to 57 km', and 84.852 to 89.716 km') consist of horizontal straight-line segments, while 
the remainder of that graph consists of 5 curved segments. The scale height, which is direct- 
ly proportional to T, is seen to remain within the range of 5.45 and 7.95 km', with a median 
value of 6.5 km'. 


*U.S. Standard Atmosphere, 1975, Committee for the Extension of the Standard Atmosphere, GPO, in preparation. 
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Figure 1. The U.S. standard atmosphere as represented by graphs of five atmospheric properties versus 
height: that is, temperature, scale height, natural log of density ratio, natural log of pressure ratio, and 

derivative of natural log of pressure ratio with respect to height, plus a sixth property for comparison- 
natural log of density ratio versus height for an isothermal atmosphere. 
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Density-Height Relationship 1 


When the differential form of the hydrostatic equation is combined with the differential 
form of the equation of state to eliminate the differential of pressure, we obtain the follow- 
ing relationship between temperature and the density-height function: 


T = 


-dh 

d£np 


"GM dT 
R + dh 

L J 


(7) 


In the absence of an analytical expression for p(h), a numerical approximation for dh/d£np 
leads to a mean value of temperature T for the height interval h 2 - hj : 


T = 


h 2 " h l 
£npj - £np 2 


GM 


T2-T, 

h 2 - hj 


In terms of scale height, this expression becomes 


( 8 ) 


H = 



Cnpj - Cnp 2 




(9) 


It is important to note that while equation (8) expresses the mean temperature for the layer 
hj to h 2 , the expression involves the average gradient of T with respect to h within that 
layer. A similar situation prevails for the expression of the mean scale height in equation 
(9) where H = (H 2 + Hj )/2. Since prior knowledge of the value of neither (T 2 - Tj )/(h 2 - hj ) 
nor (H 2 - Hj )/(h 2 — hj ) is generally available, equations (8) and (9) are not of themselves 
useful relationships for temperature-height determination. However, one version of the 
following expression derived from equation (9) will be of special importance in an uncertainty 
analysis in the section on Uncertainty of Temperatures Deduced from Densities. This 
expression is 

1 _ (H 2 +Hj)/2 

Knpj-Cn p 2 (h 2 -hj) + (H 2 -Hj) (10) 


Equation (10) is written for the case in which h 2 > hj , as when data are being analyzed from 
lower to greater heights, such that p > P 2 . If the equation is to apply to data being analyzed 
from greater to lower heights, as is the preferred case for the integral form of the density- 
height function involving the normal atmosphere, hj is greater than h 2 such that p 2 > p t , and 
equation (10) must be rewritten as 


1 (Hj +H 2 )/2 

Cnp 2 -Gnpj (hj -h 2 ) + (Hj -H 2 ) 



(10a) 


In this expression Ah is always positive, and AH is positive in regions of positive temperature 
gradients, zero in regions of zero temperature gradients, and negative in regions of negative 
temperature gradients. 
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Density-Height Relationship 2 

The integral of the hydrostatic equation, when combined with the equation of state to elim- 
inate the pressure, yields the following directly useful integral expression relating tempera- 
ture to density: 


T 


q 



1 f h q 

. — • I p(h) dh 
P q \ 


(ID 


The temperature T q is not the mean value of temperature for some layer, as in the case of 
T in equation (4), rather, it is the temperature at the specific height h q . 

It has been demonstrated (Minzner et al., 1964 and 1966) that, when applied to a helium 
atmosphere, the properties of equation (11) differ markedly from those observed when the 
equation is applied to an argon atmosphere. With an air atmosphere having a mean molecular 
weight of about 29, the properties of this equation are similar to those found when the 
equation is applied to an argon atmosphere in which the integration optimumly proceeds 
from the greatest to the lowest altitude of the density-height data. In this situation, Tj , 
pj , and hj are each associated with the greatest height of the data set, while T q , p q , and h q 
are each associated with the running value of h. This is the height for which the value of 
T q is being computed, a height which varies progressively from h t to the lowest height of 
the data set as the calculation of the profile proceeds. 


Because p(h) is usually known numerically rather than analytically, it is convenient to re- 
place the integral of equation (11) with an appropriate series approximation, one of which 
is governed by the logarithmic trapezoidal rule (Minzner et al., 1965). With the use of this 
approximation, equation (11) may be rewritten as 


P, CM 1 Ji, (hj . , -hj) (Pj -Pj . i) 

T n = — _ t, + • _ . V J - - ( 12 ; 

Pq R Pq j = 2 £n Pj - £n Pj _ j 

The logarithmic-trapezoidal rule is a particularly suitable one for the integration of atmos- 
pheric density with respect to height because this approximation very closely represents the 
conditions of the real atmosphere, and because, in a graph of measured values of £np versus 
h, the straight-line segments between successive data points exactly represent the conditions 
of the logarithmic trapezoidal rule. The agreement of this rule with reality is evident from 
the nearly straight-line segments of an idealized version of the real atmosphere shown in the 
graph of £n(p/p r ) versus height in figure 1 . For this graph, which represents the densities of 
the U.S. Standard Atmosphere, 1975, the value of p r was chosen to be I X 10' 5 kg m' 3 , so 
that the single ordinate scale at the left of the figure applies to both £n(p/p r ) and £n (p/p r ). 
The shape of the curve £n (p/p r ) versus h is identical to the shape which the curve £np versus 
h would have, because the values of £n (p/p r ) versus h and of £np versus h are offset by the 
constant difference Cnp r at all heights. 


6 



Figure 1 also contains a graph of fin (p./p r ) in the form of a continuous straight line immed- 
iately above the graph of fin (p/p r ). This single straight line diagonally across the entire figure 
is characteristic of the density of an isothermal atmosphere extending upward from a height 
of 1 1 km' where the standard-atmosphere density is 5.4721 X 10" 2 kg m' 3 , and where 
the temperature and the corresponding scale height are 216.650 K and 6.363 km', respectively. 
For this specialized atmosphere, the logarithmic-trapezoidal rule represents a series which 
exactly duplicates the integral of equation (11). The small deviations of the slopes of fin (p/p r ) 
versus height from those of fin(p./p r ) versus height show the small influence of the variation 
of atmospheric temperature from the fixed value 216.650 K (in the height region of 10 to 
90 km') on the general shape of the curves of fin (p/p r ) and finp versus h. (The influence of 
temperature-height variation upon the curves of fin(p/p r ) and finp versus h is similarly small.) 

In the evaluation of T q by equation (11), or by any appropriate approximation of that 
equation as exemplified by equation (12), a knowledge of the initial temperature T x is of 
importance only for the upper regions of the profile. For h q = h x , the value of (p x /p q ) T x 
is exactly T x because the integral term is zero and p x = p . As h q decreases from h x (that 
is, as the value of (h x - h ) increases), the value of the relative contribution of the density- 
ratio term decreases, while that of the integral term correspondingly builds up. When 
(hj - h ) increases from zero, first by one scale height and then by three scale heights, the 
relative contribution of the density-ratio term to T q decreases from 100 to about 37 percent, 
and then to about 3 percent, because of the large decrease in the value of (p x /p q ) over these 
height regions. Simultaneously, the value of the integral term grows correspondingly. At 
heights of more than three scale heights below h x , the value of T q is determined almost 
completely by the integral term alone. Figure 2 shows the relative contributions made to 
T q by the density-ratio term and by the integral term of equation (1 1) as a function of the 
range of the limits of integration (h x - h q ) expressed in units of scale height. 

The graphs in figure 2 are based on the density-height profile for an isothermal atmosphere 
(T = 216.650 K) for which finp is a linear function of height as previously shown in the graph 
of fin(p./p r ) in figure 1. Only small variations from the values of the two terms of equation 
(12) depicted in figure 2 would be seen for calculations based on a real, variable-temperature 
atmosphere. 

The elimination of T between equation (12) and equation (3) yields the following expres- 
sion for scale height: 


Hq = — * H, + • £ 


j = 2 


(hj . j - hj) (Pj - Pj - i ) 
finpj - finpj . j 


(13) 


The relative contribution of each of the two terms of this equation to H q , as a function of 
the range of integration (h x - h q ), follows exactly the same pattern shown in figure 2 with 
regard to T q . 
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ALTITUDE RANGE OF INTEGRATION (h, - h q ) (in units of one scole height) 


Figure 2. Relative contribution made to the calculated temperature by the density-ratio term and by the 
integral term of equation (11) as a function of the size of the altitude range of integration, for a 217 K 
isothermal atmosphere. 

Atmospheric soundings of both pressure and density have for many years served as the basis 
for the determination of the height profiles of temperature and scale height. These profiles 
for temperature and scale height are obtained from pressures through equations (4) and (5), 
respectively, while from densities they are obtained through one or another version of 
equations (12) and (13), respectively. It is obvious that the uncertainties of these derived 
temperatures and scale heights are a function of the uncertainty in the measured quantities 
of pressure or density. Somewhat less obvious is the fact that the height interval between 
successive observations of pressure or density strongly influences the propagation of the 
observational uncertainties into the computed temperatures and scale heights. A rigorous 
error analysis of these two pairs of equations demonstrates this situation. The pair of 
pressure-related equations are analyzed in the following section while the pair of density- 
related equations are analyzed in the subsequent section. 

UNCERTAINTY OF TEMPERATURES DEDUCED FROM PRESSURES 

The error-analysis method employed is the first-order Gaussian method, wherein each 
variable y. entering into the expression of a particular function of these variables x(y.) is 
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assumed to have an observational uncertainty 5y. which meets the conditions of a Gaussian 
or normal distribution about y { , where y. is the mean of a set of individual observations of the 
i th variable or the true value of the ith variable. Thus, if the value of x is determined from 
the functional expression x(y.), the value of 5x, the implicit uncertainty in x, is given by 

=£ fe^i ) 2 ° 4) 

Applying this relationship to the variables of equation (4), in which (Ah) t 2 is redesignated 
simply as Ah, yields 

< is > 

It may be shown from (5) that the partial derivative of T with respect to Ah, when simul- 
taneously multiplied and divided by equivalent portions of (6), leads to 


3T kAh 6 Ah 

5 Ah = • — - 

3 Ah (2np, -£np 2 ) 2 H 


(16) 


Because the coefficient of 6Ah/H in (16) is identically the coefficient of -Sp/Pj and of 
-<5p 2 /p, when the partial differentiation processes indicated in (15) with respect to both 
Pj and p 2 have been performed, it follows that 


5T 


kAh 

(£npj -£np 2 ) 2 



(17) 


where 8p l and 5p 2 are the pressure uncertainties of two consecutive pressure-height values, 
and 5 Ah is the uncertainty of the height interval Ah between the corresponding two pressure- 
height values. 

Dividing each side of equation (17) by the appropriate side of equation (4) yields the rela- 
tive uncertainty 


5T = i r/gpi y /spiV /gAhV 

f Cnpj -Cnp 2 [y Pl J \P 2 / \ H / 

It is convenient to assume that all the uncertainty of a pressure-height point is in the pressure 
and that the height increment is exact. In this case, the values of both 5p t and 5p 2 are 
correspondingly increased over their actual values, and the term involving 5Ah in equation 
(18) vanishes. It is further assumed that the relative uncertainty in measured values of p is 
the same at h x and h 2 ; that is, 5p A fp x =5p 2 /p 2 = 5p/p. These two assumptions permit 
equation (18) to be reduced to the simpler form, 


3 /2 

(18) 
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( 19 ) 


5T _ y/7 m 5p 

T Cnpj - Cnp 2 p 

From equation (6) it_is apparent that l/(2npj - £np 2 ) in equation (19) may be replaced by 
[H/(Ah)j 2 ] where H is the mean scale height for the layer h 2 - h t = (Ah) 5 2 . Thus, 
equation ( 1 9) may be rewritten as 



This equation indicates that 5T/T, the relative uncertainty in the mean temperature of the 
layer bounded by pressures p t and p 2 , is equal to the relative uncertainty in the pressure 
measurements times a coefficient 2 . This coefficient is seen to depend upon only two 
quantities, the mean temperature T of the layer h 2 - hj (T being implicit in H) and (1/Ah), 
the reciprocal of the thickness of the pressure-sampling interval (h 2 - h t ). Since the max- 
imum value of T in the earth’s atmosphere below h = 100 km' is less than twice its minimum 
value, and since Ah can, in principle, be made to vary over many orders of magnitude, it is 
apparent that Ah is the dominant factor in determining the propagation of (5 p/p) into 
5T/T). 

A graph of the coefficient Aj as a function of Ah for three specific values of T, 1 69. 1 0 K, 
241.57 K, and 3 14.04 K, is given in figure 3. The first and third of these temperatures were 
selected in part because they are close to the lowest and highest atmospheric temperatures 
normally observed at heights below 100 km, while the second is the mean of these extremes. 
In addition, these three values were specifically selected to correspond respectively to a 
particular set of three scale heights: ( 7/y/2 ), (10A/2), and (13A/2) km'. In each of these 
three cases, the value of A x 2 is unity when Ah is a particular integer multiple of one geo- 
potential kilometer, 7, 10, and 13 km', respectively. 

It is interesting to note that, for the entire range of normally observed temperatures at 
heights below about 100 km', the pressure-sampling interval Ah corresponding to unity for 
the coefficient A ; 2 varies between the limited range of 7 to 13 km'. The pressure-sampling 
height interval, however, may actually vary over several orders of magnitude depending upon 
the design of the measuring system. The value of the coefficient A t 2 could correspondingly 
vary over several orders of magnitude depending upon the choice of Ah. 

Concentrating on the median value of T, and allowing the pressure-sampling interval Ah to 
decrease first from 10 4 m' to 10 3 m' and then to 10 2 m', causes the value of A ; 2 to increase 
first from 1 to 10 and then to 100, respectively, such that the uncertainty in the value of T 
as expressed by equation (20) simultaneously increases by factors of 10 and 100 for a fixed 
uncertainty in the pressure. Conversely, as Ah is increased from 10 4 m' to 10 s m', the value 
of Aj 2 decreases from unity to 0.1 , such that the uncertainty in a related mean value of T 
correspondingly decreases by a factor of 10. 
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Figure 3. Value of the coefficient A x 2 in equation (20) as a function of pressure-sampling height interval 
for each of three atmospheric temperatures representing approximately the maximum, median, and mini- 
mum values, respectively, for heights below 100 geopotential kilometers. 
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Figure 3 shows that when mean temperatures are deduced for successive layers from pressure 
observations having a fixed uncertainty, finer height resolution in the resulting temperature- 
height profile is obtained at the expense of temperature accuracy, while increased tempera- 
ture accuracy is achieved at the expense of height resolution. Since the ultimate accuracy 
of pressures obtained from any high-altitude balloon-borne or rocket-borne pressure-sensing 
device is limited by considerations of various perturbing phenomena such as outgassing, 
boundary layer, and_shock wave, the value of 5 p/p has a practical lower bound. Thus, 

(5T/T) • (Ah/^/2*' H) which is_equal to Sp/p similarly has a practical lower bound, and for 
this minimum value of 5p/p, 5T/T is governed by the value of the ratio (Ah/H). 

Figure 3 shows the effect of variations of both T and Ah on the value of the coefficient 
Aj 2 = (\/T' H/Ah). By expressing Ah in multiples of one scale height, _the number of 
variables in equation (20)_is effectively reduced by one, and values of 5T/T can be plotted 
as a function of Ah = n • H (where n is the value along the abscissa) for any particular value 
of 8p/p. Figure 4 presents such a graph for each of nine values of pressure uncertainty. 

In figure 4, the straight-line graph for 1 percent uncertainty in p intersects the coordinate 
for 1 percent uncertainty in T at the abscissa coordinate value of 1.414. This same straight- 
line graph for 1 percent uncertainty in p intersects the ordinate values of 10, 100, and 1000 
percent at abscissa values of Ah equal to 0.1414, 0.01414, and 0.001414, respectively. 

This series of successively decreasing values of Ah, each being one tenth of the preceeding 
one, yields a geometric series of increasing values of 8T/T, each being ten times the pre- 
ceeding one. A similar situation prevails_for each of the eight other values of 5p/p, for which 
lines have been plotted. The value of 5T/T is obviously varied over orders of magnitude by 
comparable variations of the value of the ratio H/Ah. For height increments with a value 
of-v/2~times one scale height, however, the lines for each of the nine values of 5 p/p plotted 
show the relative temperature uncertainties to be equal to relative pressure uncertainties. 

In order to show simultaneously the small additional influence produced by the allowable 
range in the value of T as it appears both in 5T/T as well as intrinsically in H, the single- 
line graph for each value of 5p/p in figure 4 is expanded into a band in figure 5 where the 
pressure-sampling height is expressed in meters as in figure 3. 

The lower left-hand edge of each band corresponds to T = 169.10 K, while the upper right 
hand edge of each band corresponds to T = 314.04 K with other points across each band 
corresponding logarithmically to intermediate temperatures. From figure 5 one may estimate 
the percent uncertainty in particular atmospheric temperatures computed from pressure- 
height data measured with a specified uncertainty over particular height increments. 

Three sample applications of this graph are cited: In the first it is desired to determine 
the maximum pressure-gage uncertainty allowable to achieve a 1 percent uncertainty in a 
mean temperature of about 240 K for a layer thickness of 100 m'. We look for the inter- 
section of the 100 m' abscissa value with the 1 percent ordinate value, and find that it lies 
near the 242 K value of the 0.01 percent pressure-uncertainty band. Thus, pressures would 
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Figure 4. Percent uncertainty in the mean temperature of an atmospheric layer (for any one of nine 
uncertainties in the pressure-height data) as a function of the pressure-sampling height interval (the 
thickness of that layer) when that mean temperature is deduced from the atmospheric pressure at the upper 
and lower boundaries of that layer, and when the mean temperature is at or near 241 .57 K. 
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Figure 5. Relative uncertainty in temperatures (calculated from pressure-height data) as a function of the 
pressure-height sampling interval, and as a function of the value of the temperature, for each of nine values 
of relative uncertainty in the pressure-height observations. 



have to be measured with an uncertainty no greater than 0.01 percent to achieve the 
desired results. Since it is essentially impossible to achieve such a small uncertainty in 
any rocket or balloon measurements of pressure, this combination of height resolution and 
temperature uncertainty is essentially impossible to achieve from pressure-height data. 

In the second example, we assume a pressure-gage uncertainty of 3 percent and a mean 
temperature of 180 K, and we seek the temperature uncertainty associated with particular 
sampling-height increments. If the pressure sampling-height interval is 3 km, this ordinate 
value is seen to intersect the appropriate region of the band for 3 percent pressure un- 
certainty at a value corresponding to a temperature uncertainty on the ordinate scale of 
about 7.4 percent or about 13.3 K. A doubling of the layer thickness to 6 km would 
halve the temperature uncertainty. 

The third example involves the inverse problem of estimating the uncertainty in the 
computed pressure differential associated with an assumed isothermal layer of fixed 
thickness. If the layer has a thickness of 3 km, as may be the situation in the grenade 
experiment (Nordberg and Smith, 1964), and if there is an uncertainty of about 2 per- 
cent or 5 K in an assumed mean temperature of 250 K, the graph shows that the com- 
bined uncertainty of the two boundary pressures lies between 0.3 and 1.0 percent or 
about 0.6 percent on a logarithmic scale. For a layer thickness of 1 km, the boundary- 
pressure uncertainties would increase by a factor of 3. If the computed pressure-height 
profiles were used to generate density-height profiles, the uncertainties in the densities 
would be somewhat larger than in the pressures because of the required vector addition 
of the temperature uncertainty and the pressure uncertainty to obtain the density 
uncertainty. 

UNCERTAINTY OF TEMPERATURES DEDUCED FROM DENSITIES 

Applying the first-order Gaussian method, equation (14), to the determination of the 
uncertainty in temperatures deduced from density-height data through equation (11) 
leads to 



provided that we assume all the uncertainty in a density-height point to be concentrated 
in the density. Because scale height combines the temperature with several constants, 
the use of scale height facilitates the expansion and analysis of equation (21). Consequent- 
ly, it is convenient to rewrite that equation as 
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Both equations (21) and (22) involve (q + 1) terms, where q is the number of density- 
height data points. The first term is associated with the uncertainty of the temperature 
or the scale height at hj , while the second term is associated with the uncertainty of the 
density of h x . Each of the q-1 additional term deals with the uncertainty of one of the 
successive q-1 density-height data points, respectively. Three of the (q + 1) terms of the 
series have unique formats. These are the first, the second, and the last terms, those 
involving 5Tj (or SHj ), 8p 1 , and 5 p q , respectively. However, the remaining terms 
(the third through the qth term involving 5 p 2 through 5p ql , respectively), have a com- 
mon format. Thus, the sum of these common-format terms may be expressed as the sum- 
mation of a general term. Applying the operators indicated in equation (22), dividing 
both sides of the resulting equation by H q , and introducing equation (10a) yields the 
following equation: 
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Equation (23) (developed in appendix A) is very much more complicated than the com- 
parable uncertainty expression in terms of pressure-height data as represented by equation 
(18). The considerable difference in the complexity between the two equations stems 
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largely from the fact that equation (18) involves the uncertainty of only two pressure- 
height data points and the uncertainty of the corresponding height increment, while 
equation (23) involves the uncertainty of all consecutive density-height data points from 
hj to h q , each of which contributes at least a small amount to the value and uncertainty 
of the temperature or scale height at height h q . It is emphasized that the quantities T q 
or H q in equations (11) through (13), and again in equations (21) through (23), represent 
values for specific heights rather than mean values for specific layers as in the case of the 
pressure-related equations. 

Each element of equations (23) through (27) is associated either with a specific height 
or with a specific layer. Elements having a single subscript, for example, H x , p x , p., p. x , 
p j+1 , p q x , p q , and H q , signify the value of the particular quantity at heights h x , h., h jx , 
h. + x , h j , and h q , respectively. Quantities with a double subscript, for example, 

«i 2 , ce. j ., and so forth, represent a quantity associated with particular layers, that is, 
the layers bounded by h x and h 2 , or by h. x and h., and so on, respectively. Thus a x 2 , 
as expressed by equation (24),_represents an algebraic expression of four different quan- 
tities, three of which, namely H, Ah, and AH, are associated with the layer h x to h 2 . 

In equation (24), H is the mean scale height for the layer Ah = h x - h 2 , while AH 
represents the change in scale height within that layer, that is, H x - H 2 , such that only 
in a region where the gradient of H or T with respect to height is positive will AH be 
positive. The quantity H q is, of course, the value of H at height h q . Similarly, a. x 
a. j+1 , and a , as defined by equations (25), (26), and (27), respectively, each 
represents quantities associated with the particular appropriate layers. Equation (23) 
also includes the doubly subscripted quantities 

( J._) (_2_) ( H ( R ) 

\Ah + AH/ j 2 \Ah + AH/j _ j j 5 \Ah + AH/^ ■ + , ' \Ah + AH/ q _ x ? q 

each of which represents one factor of the right-hand side of equations (24), (25), (26), and 
(27), respectively. The influence of nonzero gradients in the temperature-height profile 
is impressed upon 5H q /H q through all the doubly subscripted quantities in equation (23). 

In order to see more clearly the influence on 5H q of sampling-height interval alone, it is 
convenient to assume an isothermal atmosphere and to correspondingly simplify equation 
(23). In this case, scale height would not vary but would remain fixed at a value H over 
the_entire region of integration, and the following relationships would apply: H 1 = H q = 
H, H = H, and AH is zero in all layers. Under these conditions, equations (24), (25), 

(26), and (27) become unity, as does the ratio Hj/^, and equation (23) may be rewritten 
as 
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It can be shown that 



(29) 


such that these two terms in the coefficient of (8p q /p q ) in equation (28) cancel each 
other. Then, if we impose the additional condition that the height increments between 
successive density-height data points are constant over the entire region of integration, 
this condition plus equation (29) permit the further simplification of equation (28) to 
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(30) 


If we impose still another restriction, that is, that the relative uncertainty of the density 
data has the constant value 6p/p for all heights within the range of integration, we see 
that Spj/Pj = 5 p./p. = 5p q /p q = 8p/p. Thus we may rewrite equation (30) as 



The restriction permitting this simplification, while somewhat unrealistic from the point 
of view of any measuring system, is acceptable because the relative uncertainty of only 
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p q enters significantly into the value of 8H q /H q as shown below. It is convenient to 
represent the sum of the series of the squares of the several density ratios in the coefficient 
of (8p/p) 2 in equation (31) by (S q ) 2 , such that 



With this simplification, equation (31) may be rewritten as 



(33) 


Finally, recalling from figure 2 that the ratio (P x /P q ) causes the contribution of to H q 
in equation (13) to become negligible when the integration has proceeded downward from 
hj by more than about three scale heights, we recognize that the contribution of SHj/Hj 
to 8H q /H q in equation (33) must similarly be negligible for h q sufficiently below hj . 

Thus for this condition, equation (33) may be approximated by 



(34) 


Also, because 8H q /H q is identically equal to 8T q /T q , it is convenient at this point to 
return from scale-height notation to temperature notation, and equation (34)* is rewritten 
as 



(35) 


This equation represents the relative uncertainty in T q (the temperature at height h q ) as 
deduced from a set of density-height data measured with a constant relative uncertainty 
(8p/p) over the entire range of a height region h t down to h q , where (h x - h q ) is equal 
to at least three scale heights. _This equation assumes the atmosphere to be isothermal 
at the value T associated with H, a condition which might apply to the earth’s atmos- 
phere at heights above about 400 km. Because several of the variables included in the 
most general form of the uncertainty expression as given by equation (23) are excluded 
in equation (35), this simpler version is more readily examined for the influence upon 
8T q /T produced by variations in the sampling-height interval Ah alone. 


*At this point and in the balance of the paper the word equation shall be construed to include mathematical statements 
of approximate equalities such as in (34). 
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In terms of the relative_uncertainty in density, equation (35) for 5T q /T q is very similar 
to equation (20) for 5T/T in terms of the relative uncertainty in pressure, except that S q 
in equation (35) takes the place of ^/Tin equation (20). Therein lies a_great difference, 
because the quantity S q is a dimensionless function of three variables, H, Ah, and 
(hj - h q ), and has a value which varies between y/T and zero, as Ah varies from large 
values to zero. For Ah > 5 H, the value of S q approaches y/T asymptotically as Ah 
approaches infinity, independent of the value of H or (hj - h q ). For Ah < 0.2 H, the 
value of S decreases nearly linearly with decreasing values of Ah, for any fixed value of 
H or (hj -h q ). 

The influence of variations of both Ah and (h x - h q ) upon S q , and implicitly upon 
5T q /T q , is seen in figure 6. Here, for a fixed value of H = (10 /y/T) km, a band of values 
of S q is plotted as a function of Ah on a fully logarithmic scale. Ordinate values of S q 
are given on the right-hand side of the figure. The upper limit of the band is designated 
S q max and is plotted as a line of long dashes; the lower limit of the band is designated 
S q ^ and is plotted as a solid line. The range of the values of S q within the band, at 
any particular value of Ah, represents the influence of variations of (h t - h q ) at that 
value of Ah, with the value of S q increasing as (h t - h q ) decreases. At any particular 
value of Ah, the smallest possible value of (h t - h q ) is, of course, Ah, and hence S q max 
depicts the locus of the values of S q associated with the condition that (h t - h ) = Ah. 
Thus, this upper limit of S q represents the value of equation (32) when that equation 
involves only the first two density-altitude points of a data set, that is, the value of S q 
for the case when 5T q /T q is being computed for the height h q associated with the 
second-highest density-height data point in the sounding. For this situation, the right- 
hand side of equation (32) reduces to [y/2 (p q - Pj)/p q ]. As the height for computing 
5T q /T is lowered from hj , such that (hj - h q ) becomes increasingly large, the value of 
S q slowly decreases and asymptotically approaches a lower limit S q min when (h t - h q ) 
is greater than three scale heights. 

It is convenient to define the ratio of S„ to S„ m . as R and to examine the variation 
of R with respect to variations in Ah. As long as Ah is smaller than about 0.2 H, R 
retains a nearly constant value (a value which is seen graphically to be approximately 
equal to y/T), and the graphs of both S„ and S„ . versus Ah are seen to make an 
angle of about +45° with respect to the abscissa. Because of the feature of a nearly 
constant slope of +1 for both S„ and . in the region where Ah < 0.2 H, this 
region is hereafter designated as the small-increment regime. As Ah exceeds 0.2 H and 
increases toward 5 H, the slopes of both S„ and S„ m . gradually decrease toward 
zero. However, the slope of S decreases more rapidly than that of S .so that 

5 x q,max r J q, mm 

S approaches S . at a common value yfl, while the value of the ratio R decreases 

q, max r q, mm v 

from s/2 toward unity. For Ah > 5 H, the value of R is essentially constant at unity, 
and the slope of the line common to both S„ and S„ . is essentially zero. The 
feature of an essentially zero slope for this common line representing S q versus Ah specifies 
a region which, because of the related large value of Ah (Ah > 5 H), is hereafter referred 
to as the large-increment regime. The region between the small-increment regime and the 1 
large-increment regime, the region for which 0.2 H < Ah'< 5 H, is hereafter referred to as i 
the transition regime. 
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Figure 6. Graphs of three quantities as a function of density-sampling-height interval Ah, that is, (1) the 
limited band of values of a particular series of terms defined by S q , (2) the value of the last term of that 
series, and (3) the value of the ratio of mean scale height to Ah. 

In the large-increment regime, where S q max and S q min are essentially identical, it is 
apparent that S q is essentially independent of (h 1 - h q ). In the small-increment regime, 
the concern for the influence of (h 1 - h q ) upon S q remains only for values of (h x - h q ) less 
than about 3 H. Because equation (35) is based upon the restriction that ( h l - h ) is greater 
than 3 H, and because S q approaches S q min for such values of ( h 1 - h q ) in all three regimes, 
it is immediately apparent that the general factor S q in equation (35) should be replaced by 
the specific value S q min . It is somewhat less apparent that 5p/p (the general expression for 
relative uncertainty of density in equation (35)) should be replaced by 5p q /p q (the specific 
relative uncertainty for h q ), particularly when equation (35) is applied to the small-increment 
regime. The reasons for the latter replacement stem from the facts that the following dual 
relationship exists for the dual conditions (h 1 - h q ) > 3 H and Ah < 0.2 H: 

• The value of the series S q , which becomes essentially S q min for ( h 1 - h q ) > 3 H, 

simultaneously approaches the value of (p q -p )/p , the last term in equation (32), 
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defining S q . This near equality between S q mjn and (p q - p ql )/p q in the small- 
increment regime is demonstrated graphically in figure 6. 

• The term (p q - p q l )/p q is properly associated with the specific uncertainty 
5p q / p q rather than with 5p/p. This situation is evident from the fact that the 
quantity (p q - p q l )/p q is a coefficient of 5p q /p q in equation (30), which, 
except for an isothermal condition, is a general expression equally applicable to 
both 8H q /H q and 5T q /T q . Thus, at least for the small-increment regime, equation 
(35) may be rewritten as 



(36) 


Even in the large-increment regime where the value of S to approaches [\Z?(p q - P q . 1 )/p q ] , 
the value of S q mjn is still dominated by (p q - p q l )/p q (which is equal to unity in this 
regime); therefore, S . should still be associated with 8p n /p . 

0/7 q, min q q 

Because (p -p x )/p approaches unity in the large-increment regime, while S q min ap- 
proaches •s/7, it is reasonable to write the following special form of equation (36) for that 
regime : 



This equation is seen to be analogous to equation (20) involving uncertainties in pressure- 
height data. 

In addition to the graphs of S q>max , S q mjn , and (p q - p )/p q versus Ah, figure 6 also 
contains a graph of the function (H/Ah) versus Ah for H = (10/yT) km or for T = 241.57 K. 
This function, which obviously varies inversely with Ah, is seen to have a constant negative 
45° slope for a fixed value of H. When this function is multiplied by S q min , a quantity 
which according to figure 6 is seen to have a +45° slope in the small-increment regime, the 
product for the small-increment regime is essentially a constant (independent of Ah) with 
a value near unity for all values_of H or T. This situation is depicted in figure 7 where the 
value of the product [S q min • (H/Ah)] versus Ah is plotted in all three regimes for each of 
three values of H — (7/^/7), (10A/2), and (13/y/J) km' — consistent with the mean tempera- 
tures 169.10, 241.57, and 314.04 K, respectively. Because of the characteristics of this 
product, equation (35) as applied to the small-increment regime may be replaced by 



( 38 ) 
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Figure 7. Value of the product of two factors (comprising the coefficient of the uncertainty of observed 
densities) in equation (36) as a function of density-sampling height interval, for each of three values of 
mean scale height consistent with three specified temperatures. 

Thus, three different expressions represent the value of 5T q /T q as a function of uncertain- 
ty in density-height data: equation (38) in the small-increment regime, equation (37) 
in the large-increment regime, and equation (36) in the transition regime. 

As a representative of the three equations noted, the application of the data from figure 7 
to each of nine particular values of uncertainty in density-height data yielded the graph 
shown in figure 8. This figure shows the percent uncertainty in the temperature, derived 
from density-height data through equation (12), as a function of density-sampling height 
interval for a band of normal atmospheric temperatures, 169 to 314 K, and for a wide range 
of density-height uncertainty, 0.01 to 100 percent. This figure can serve as the basis for 
estimating the relationship between the uncertainty in each value of a set of density-height 
data to each point of the related temperature-height profile, or vice versa, for the region of 
the earth’s atmosphere below about 120 km. 

The examples of hypothetical use of the pressure-data graph, discussed in the section 
on temperatures deduced from pressures, could be readily transposed to apply to this 
density-data graph. In particular, figure 8 shows that the relative uncertainty in the 
density observations need not be smaller than 1 percent in order to have less than a 1 per- 
cent uncertainty in a 245 K temperature for a temperature height profile with a height 
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Figure 8. Relative uncertainty in temperatures (calculated from density-height data) as a function of the 
density-sampling height interval and as a function of the value of the temperature for each of nine values of 
relative uncertainty in the density-height observations. 



resolution of 100 meters. This value is 100 times larger than the 0.01 percent uncertainty 
required of pressure data to meet the same conditions and represents a very significant re- 
laxation of the measurement requirements. 

The estimates deduced from the use of this graph are theoretically accurate as long as the 
data to which the graph is applied comply with the three restrictions imposed during the 
preceding development. However, it will be seen that one of these restrictions may be 
essentially eliminated, while a second has only a small influence. One may recall that these 
restrictions are: 

• The height h q associated with T q is more than three scale heights below h x the 
greatest height of the sounding: (h x - h q ) > 3 H, 

• The relative uncertainty of the density-height data is constant over the entire 
height range of the sounding, and 

• The atmosphere is isothermal over the height range of the sounding. 

The first of these restrictions must be retained unless one has some independent means for 
determining Tj or Hj and its uncertainty as required by equation (33). Only for (hj - h )> 
3 H will the ratio Pj /p q be sufficiently small to permit the associated term to be neglected 
in equations (34) through (38). 

The second restriction is not a very significant one because it has been shown that the 
constant factor (5p/p), used as the general expression for density uncertainty in equation 
(35), is dominated by the particular value (5p q /p q ) associated with the density-height data 
at h q . Thus, the need for (5 p/p) to be constant over the entire height range of a particular 
sounding can be relaxed without significantly affecting the validity of any of equations (36), 
(37), or (38). 

The third restriction appears to be philosophically important because the earth’s atmosphere 
below 120 km is certainly not isothermal. Actually, however, the existence of the various 
nonzero gradients in the temperature profile has little effect on 5T q /T q . This situation is 
due to the fact that 5T q /T q depends primarily on the conditions between the data points 
(p q t , h q j ) and (p q , h q ) and hardly at all upon the conditions between other pairs of data 
points. Even a nonzero temperature gradient between the two specified data points has 
only a small effect, and this can be readily accounted for. This is accomplished by reintro- 
ducing into equation (36) certain temperature-gradient-dependent coefficients which are 
associated with 5p q /p q in the general version of 5H q /H q (or equivalently 5T q /T q ), as ex- 
pressed by equation (23). 

To begin with, (H/Ah) in equation (36) is replaced by its original form H/(Ah + AH) as used 
in equation (23). Then S q min , which has been shown to be almost exactly equal to [(p q - 
Pq-i )/P q 1 f° r isothermal conditions in the small-increment regime, is replaced by the prod- 
uct of this ratio times its associated coefficient a t from equation (23). Thus, for an 
atmosphere with varying temperature or scale height, we may rewrite equation (38) for 
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the small-increment regime as 



Using the value of a q x as defined by equation (27) and replacing H by (2H q + AH)/2 where 
AH = (H -H q ), and remembering that Ah = (h q t -h q ), we have an expression for 
5T q /T q in terms of H q x and p ql at height h q t , and in terms of H q , p q , and 5p q at height 



It is interesting to examine an example of the use of equation (40) in terms of a realistic 
set of data, applicable to the small-increment regime, as taken from the U.S. Standard 
Atmosphere, 1962. Choosing the height interval Ah = 1 km' between h q x =66 km' and 
h =65 km', in a region of negative temperature gradient AT/Ah = -4 K/km', we find p , 
and p q to be 1.3482 X 10' 4 and 1.5331 X 10‘ 4 kg m" 3 respectively, H q = 7.0709 km', and 
AH = -0.1 173 km'. The substitution of these data into equation (40) leads to a value of 
1.076 for the coefficient of 5p q /p q , a value involving only a second-order difference from 
the unity coefficient of equation (38); in this case, 5T q /T would be 1.076 times 5p q /p 
instead of 1 .000 times 5p q /p for the isothermal case. Thus the assumption of isothermality, 
both in the development and use of equations (36), (37), and (38), is seen to introduce only a 
second-order error, and accounting for this error is unnecessary in most uncertainty deter- 
minations. 


CONCLUSIONS: A COMPARISON OF THE DENSITY-DATA NOMOGRAM WITH THAT 
FOR PRESSURE DATA 

The graph relating 5T/T to sampling-height interval of density data for various temperatures 
and density uncertainties depicted in figure 8 is based on the same range of values of tempera- 
ture and sampling-height interval used in figure 5, which relates 6T/T to the sampling-height 
interval of pressure data. In addition, the nine assumed values of 5 p/p in figure 8 are iden- 
tical to the nine assumed values of 5p/p in figure 5. Thus, the two figures are exactly 
comparable, and the only difference is that figure 8 is for density data while figure 5 is for 
pressure data. A comparison of these two figures shows that, at least from the point of view 
of height resolution and uncertainty of derived temperatures, density-height data are to be 
preferred over pressure-height data. 

In the large-increment regime, Ah > 5H, figure 8 is essentially identical to that of figure 5. In 
this regime, both figures show that the uncertainty in the derived temperature is decreased 
as the height resolution of the related temperature profile is decreased, that is, as the 
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sampling-height interval is increased. However, even the smallest sampling-height interval 
within this regime is so large (Ah = 5H or about 30 km'), that the height resolution makes 
this regime essentially useless for temperature-height determinations at heights below 
100 km'. 

In the small-increment regime as well as in the transition regime, figures 5 and 8 are quite 
different from each other. Assuming that (h t - h q ) > 3H, a basic condition for the appli- 
cability of figure 8, the uncertainty in temperatures derived from density-height data is seen 
to be independent of both Ah and T in the small-increment regime (Ah < 0.2H). In this 
regime, 5T q /T q derived from densities is dependent only upon the uncertainty of the den- 
sity-height value at height h q . It is apparent that uncertainty considerations place no limits 
on the usable fineness of the height resolution of the density-height data. This situation is in 
contrast to that associated with pressure-height data where ST /T is seen in figure 5 to depend 
upon T and Ah, as well as upon Sp/p, and where 5T /T becomes prohibitively large for 
reasonable values of Sp/p when Ah becomes smaller than about 1 km. 

In the transition regime, figure 8 shows less difference from figure 5 than in the small-incre- 
ment regime. In this transition regime, the characteristics of the function determining 
ST/T from density-height data vary between those of the small-increment regime and those _ 
of the largedncrement regime, such that ST/T decreases slightly as Ah varies from about 0.2H 
to about 5H. In going from the small-increment regime to the transition regime, however, 
the increased coarseness of the height resolution of the related temperature-height profile 
would more than offset the correspondingly small decrease in temperature uncertainty, 
particularly since the minimum values of Ah in this regime are already of the order of 1 km'. 
Even in this regime, however, the density-height data yield smaller values of ST/T than are 
obtained from pressure-height data. 

In general, a comparison of figure 8 with figure 5 shows that the density-height data are 
far more desirable than pressure-height data at least from the point of view of the size of 
the uncertainty and of the height resolution of the derived temperature-height profile. 


Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland August 1975 
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APPENDIX A 

DEVELOPMENT OF THE GENERAL EXPRESSION FOR 
(5 T q /T q ) 2 AND FOR (5 H q /H q ) 2 

It is well known that temperature may be deduced from density-height data by means of the 
following relationship: 


T 


q 



GM 

IT 



p (h) dh 


(A-l) 


where 

p(h) is the function relating atmospheric density to geopotential height, 

hj is the geopotential height of the upper limit of the region of integration in geo- 
potential meters (m'), 

h q is the geopotential height of the lower limit of the region of integration in geopotential 
meters (m'), 

G is the geopotential gravity constant 9.80665 m 2 s' 2 (m')' 1 , 

M is the mean molecular weight of the air 28.9644 kg (kmol)' 1 , 

R is the universal gas constant 8.31432 x 10 3 joules K' 1 kmol" 1 , 

Pj is the atmospheric density at hj , 

Tj is the atmospheric temperature at hj 

p q is the atmospheric density at h q , and 

T q is the sought after atmospheric temperature at h q . 

Because of the defined relationship between temperature and scale height H, H = TR/GM, 
it is convenient to rewrite equation (A-l) as 


Pi 1 

'’< h)dh 


(A-2) 


where 

Hj is the scale height at hj , and 
H q is the scale height at h q . 
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Generally, p(h) is not known as an analytical function for which one might find a perfect 
integral, but rather is known as a set of numerical values of density versus geometric or 
geopotential height. Therefore, it is convenient to replace the integral term in equation 
(A-2) with a series approximation. One possible form of such an approximation is a series 
of terms each of which represents the area of a trapezoid corresponding to the area under 
that portion of the graph of the natural logarithm of density versus h represented by two 
successive density-altitude points. In such a situation, the sum of these terms represents the 
sum of the areas of all of the successive trapezoids between the specified limits. When the 
successive density-altitude points plotted on a semilogarithmic scale are connected by straight- 
line segments, as is frequently the situation for closely spaced density-height data, the series 
of logarithmic trapezoids exactly fits the area under the graph. When using this logarithmic 
trapezoidal approximation, (A-2) may be rewritten as 


H 


q 


. H + — • Y' ( h j-i h j)(Pj ~Pj-i) 
P q 1 P q CnP j _CnP j-l 


(A-3) 


The validity of using (A-3) as an approximation of (A-2) for the case of a continuous atmo- 
sphere improves as the height increment between successive density-height values decreases. 

The uncertainty in the computed value of H q is based on a function involving the partial 
derivative of H q with respect to each of the independent variables. These include H f and 
the appropriate number of density-height data pairs: hj , p l ; hp p. (for j = 2 to q - 1); as 
well as h q , p q . With the assumption that the uncertainty in each data pair is entirely in the 
density value, we can assume fr j and h. to have no variability. Consequently we are interested 
in the partial derivatives of H q , as expressed by equation (A-3), with respect to only the 
following variables: Hj , p 1 , p. (for j = 2 to q-1) , and p q . The general member of this set of 
variables is given the general designation y.. The partial derivative 3H q /9y. multiplied by the 
corresponding uncertainty 5y. for each of the q + 1 variables is discussed below: 


The product of 5Hj times the partial derivative of H q with respect to is .simply 



5H, 



(A-4) 


The product of 8p r times the partial derivative of H q with respect to p l is 



Ah 


v Cnp 2 - £n p u 



,£n P 2 - 2np ly 


£ Pi 

p i 


(A-5) 
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It can be shown, however, that 


1 

Cnp 2 - £npj 

where the double subscript “I, 2” on the right-hand member of equation (A-6) indicates 
that Ah = (hj - h 2 ), AH = (Hj - H 2 ), and H = (Hj + H 2 )/2. When h t is greater than h 2 , 
the left-hand side of equation (A-6) is positive, Ah is positive, and AH has the sign of the 
temperature gradient 3T/3h in the region hj to h 2 . The same equation, with other con- 
secutive digits such as 2, 3, or 3, 4, and so on, indicates these same relationships for the 
corresponding height intervals. 

Equation (A-6) combined with equation (A-5) yields 


= ( " ) 
\Ah + AH ) j 2 


(A-6) 


9H q . (Pi H 1 17 H Ah \ ir Pl (P 2 -Px\( H \ "H *Pl 
9p, Pl "\p q 1 ~l\l Ah + AH; it2 J|_p q -^ p q ) \Ah + AH^ < 2 J/ p, 


(A-7) 


When modified by the introduction of equation (A-6), with each of two different but 
appropriate pairs of subscripted digits, the product of Sp 2 times the derivative of H q with 
respect to p 2 becomes 


3H 


q 


3 P 2 


5p 2 



In this equation the subscript “1, 2” on each of two factors has the same significance as it 
has on these same two factors in equationJA-7), while the subscript “2, 3” on each of two 
other factors signifies that the quantities H, AH, and Ah within each of these factors are 
associated withjthe geopotential height increment h 2 to h 3 . Thus, for these two factors, 

Ah = h 2 - h 3 , H = (H 2 + H 3 )/2, and AH = H 2 - H 3 . Equation (A-8) involves the two height 
increments which are separated by the height h 2 . 

The partial derivative of H q with respect to each of p 3 , p 4 , . . . , p 2 , and p q } , when 
multiplied by 5p. (where j is successively 3, 4, . . . , q - 2, and q - 1) and, when further 
modified by the appropriate introduction of equation (A-6), is identical to equation (A-8) 
except for the successive incrementing of the subscripts. Thus, a convenient form of the 
product of 5p times the particular expression for the partial derivative of H with respect 

to Pq.i is 
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' Ai 'Ah+A Hy / q _ lq j [ Pq \ Pq ;vAh + AH; q _ lq j/ Pq _, 


(A-9) 


In equation^ A-9) the subscript “q - 2, q - 1” on each of two factors signifies that the 
quantities H, AH, and Ah within each of those two factors are associated with the geo- 
potential-height increment h q 2 to h q t , so that Ah = h q 2 - h q x , H = (H q 2 + H q j )/2, and 
AH = H q 2 - H q j . Similarly in equation_£A-9), the subscript “q-1, q” on each of two 
other factors signifies that the quantities H, AH, and Ah within each of these two factors 
are associated with the geopotential-height increment h to h q , so that Ah = h q t - h q , 

H = (H q j + H q )/2, and AH = H q t - H q . This expression involves the two height 
increments separated by the height h t . 

Because of the common form of the partial derivatives of H q with respect to p 2 ,p 3 , , 

p 2 , and p q j , it is desirable to write a general version of equations (A-8) and (A-9) to 
express the product of 5p. times the partial derivative of H q with respect to p. , where j 
is understood to have values ranging from 2 to q - 1 . 

This general equation is 
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(A- 10) 


Finally, the product of 5p q times the partial derivative of H q with respect to p q when 
modified by the introduction of equation (A-6), is 
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In this equation the subscript “q - 1, q” on each of two factors has the same significance 
as in equation (A-9), while the subscript “j - 1, j” on another factor in the general term sig- 
nifies that the quantities H, AH, andAh within that term are associated with the general 
geopotential height increment h. , to h.. Thus, for this factor, Ah = h, - h., H = (H. + 

H.)/2,and AH = H._ 1 -H.. 

It is evident that one part of this equation deals with the data associated with the single 
height increment Ah = (h ql - h q ) between the lowest two density -height values involved 
in the calculation of H q , while another part of the equation deals with all the height incre- 
ments between h t and h q and their associated density data. 

Each of equations (A-4), (A-7), (A-10), and (A-l 1) is directly involved in the Gaussian 
expression for relative uncertainty 5H q /H q which follows: 



The first, second, and fourth terms within the brackets of this equation are seen to represent 
exactly the squares of equation (A-4), (A-7), and (A-l 1), respectively. The third term 
represents the square of equation (A-10) evaluated for j ranging from 2 to q - 1. The 
bracketed portion of the right-hand side of equation (A-12) is seen to be multiplied by the 
reciprocal of H q , thereby implying that each term of equations (A-4), (A-7), (A-10), and 
(A-l 1) must ultimately be multiplied by 1/H q , either before or after these equations are 
squared and summed to equal (6H q /H q ) 2 . It is convenient in this case to do the multiplica- 
tion before squaring and summing. Equations (A-7), (A-10), and (A-l 1), each contain the 
doubly subscripted factor 


(H * Ah \ 

\ 1 * Ah + AH 7 

in at least one term. It is convenient therefore to accomplish this multiplication operation 
in the appropriate terms by introducing H q into the denominator of this doubly subscripted 
factor, thereby converting this factor into a nondimensional coefficient with a value close 
to unity. In order to put the resulting coefficient into proper perspective with respect to 
the remainder of the equation and also to conserve space in the uncertainty expression 
being developed, it is convenient to define each of these modified factors by a specific 
symbol. Thus, from equation (A-7) the modified factor is defined as the coefficient a t , 
that is, 


JJL . 

1,2 \H q Ah + AH ) l2 


(A-l 3) 
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while the two modified factors involved in equation (A-l 1) are defined as 


Vij 



Ah 


• Ah + AH 


[ )hj 


(A-l 4) 


and 


a- 


jj+l 


/ _H Ah \ 

\ q * Ah + AH^j j+1 


(A-l 5) 


the modified factor from equation (A-10) is defined as the coefficient a , that is, 


a. 


q-i.q 


[ H Ah \ 

\Hq * Ah + AHy q _j q 


(A-l 6) 


Substituting equations (A-4), (A-7), (A-10), and (A-l 1), respectively, into the successive 
terms on the right-hand side of equation (A-l 2), dividing each of these equations by H q 
(that is, replacing each of the “l’s” in the denominators of these equations by H q ), and 
simultaneously replacing the resulting modified doubly subscripted factors in these equations 
by the equivalent coefficient forms defined in equations (A-l 3) through (A-l 6) leads to the 
following expression for (5H q /H q ) 2 : 



Because 5H q /H is identically equal to 5T q /T q , equation (A-l 7) may be rewritten as 
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